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We report on muon spin rotation (/xSR) measurements of the internal magnetic field distribution 
n{B) in the vortex solid phase of YBa2Cu30;; (YBCO) single crystals, from which we have simul- 
taneously determined the hole doping dependences of the in-plane Ginzburg-Landau (GL) length 
scales in the underdoped regime. We find that Tc has a sublinear dependence on 1/A^f,, where Xab 
is the in- plane magnetic penetration depth in the extrapolated limits T^O and H^O. The power 
coefficient of the sublinear dependence is close to that determined in severely underdoped YBCO 
thin films, indicating that the same relationship between Tc and the superfluid density is maintained 
throughout the underdoped regime. The GL coherence length ^ab (vortex core size) is found to in- 
crease with decreasing hole doping concentration, and exhibit a field dependence that is explained 
by proximity-induced superconductivity on the CuO chains. Both Xat and ^at axe enhanced near 
f/8 hole doping, supporting the belief by some that stripe correlations are a universal property of 
high- Tc cuprates. 



PACS numbers: 74.25.Qt,74.72.Bk,76.75.-f i 
I. INTRODUCTION 

Abrikosov vortices in a superconductor are governed 
by two characteristic length scales. The core of a vor- 
tex has a size dependent on the Ginzburg-Landau (GL) 
coherence length ^, while the supercurrents circulating 
around the core decay on the scale of the GL penetration 
depth A. In the early days of high-Tc superconductiv- 
ity it was common practice to infer the behavior of the 
in-plane magnetic penetration depth Xat from measure- 
ments of the muon depolarization rate a in the vortex 
state of polycrystalline samples. The temperature de- 
pendence of a was found to be consistent with s-wave 
pairing symmetrjsii^ii^ii^ and a universal linear scaling of 
Tc with a was observed in the underdoped regime (the 
so-called 'Uemura plot'), indicating that Tc (xl/A^^cxps, 
where ps is the superfluid density.AiS Later, microwave^ 
and ^SR^ measurements on YBa2Cu30y (YBCO) single 
crystals in the Meissner and vortex phases established 
a limiting low-temperature linear T dependence of Xat 
that is consistent with c?-wave pairing. More recently, 
non-/iSR studies of YBCO in the Meissner phase have re- 
vealed that Tc has a sublinear dependence on l/A^fci^ d^d^ 
On the other hand, the relation Tc oc ps inferred from the 
Uemura plot is supported by a recent study of electric- 
field induced superfluid density modulations in a single 
underdoped ultra-thin film of La2-a;Sr2;Cu04>^ 

The problem with assuming a oc l/A^^ is that there 
are additional inseparable contributions to a from elec- 
tronic magnetic moments and flux-line lattice (FLL) dis- 
order, which may vary with doping. To circumvent this 
difficulty we have studied YBCO single crystals. In a sin- 
gle crystal the FLL contribution to the /iSR line shape 



n{B) is asymmetric and distinct from the other sources 
of field inhomogeneity.^'^ Not only can the behavior of 
Xab be isolated, but because the finite size of the vortex 
cores is apparent in a single-crystal measurement of n{B), 
S,ab can be simultaneously determined. While ^ab may 
be accurately determined in conventional superconduc- 
tors from measurements of the upper critical field Hc2, 
in high-Tc cuprates Hc2 is generally a very high mag- 
netic field marking the transition from a vortex liquid to 
the normal phase. Here we present /iSR measurements 
that probe Xab and £_ab in the bulk of YBCO single crys- 
tals deep in the superconducting state. The accuracy of 
our method was demonstrated in previous studies of con- 
ventional superconductorsr^ii^ii^ but is reinforced here 
through comparisons with the results from other tech- 
niques. 



II. EXPERIMENTAL DETAILS 

YBCO single crystals with purity of 99.995 % were 
grown by a self-flux method in fabricated BaZrOa cru- 
cibles at the University of British Columbiaj^i An ex- 
ception are y = 6.60 single crystals that were grown in 
yttria-stabilized-zirconia crucibles and characterized by a 
purity greater than 99.5 %.— Single crystals of Ca-doped 
YBCO were also prepared in BaZrOa crucibles. Typical 
sample sizes consisted of 3 to 5 single crystals from the 
same growth batch arranged in a mosaic to form a total 
a-b surface area of 20-30 mm^ . The thickness of the crys- 
tals are on the order of ~0.1 mm. The superconducting 
transition temperatures of the single crystals were mea- 
sured using a SQUID magnetometer. Twin boundaries 



TABLE I: Characteristics of the YBa2Cu30y and (Y, 
Ca)Ba2Cu3 06.98 single crystals. The hole concentration p 
per Cu02 layer is determined from the dependence of Tc on 
p presented in Ref.— for similar YBCO single crystals. 



y 



p Tc (K) Detwinned 



6.60 
6.57 
6.67 
6.75 
6.80 
6.95 



0.103 
0.110 
0.120 
0.132 
0.140 
0.172 



62.5 
59.0 
66.0 
74.6 
84.5 
93.2 



V 
V 
V 



(Y,Ca)6.98 0.192 86.0 



were removed from some of the higher-doped single crys- 
tals by applying pressure along the a or 6 directions at 
elevated temperature. These basic sample characteristics 
are summarized in Table IT] 

The ^SR experiments were performed over a 3 year 
period on the M15 and M20B surface muon beam Hnes 
at TRIUMF, Vancouver, Canada. The /iSR spectra were 
recorded in a transverse-field (TF) geometry with the ap- 
phed magnetic field H perpendicular to the initial muon- 
spin polarization direction, and perpendicular to the c- 
axis of the single crystals. A TF-/iSR spectrum com- 
prised of 20 to 30 million muon decay events was taken 
at each temperature and magnetic field. 

In a trans verse- field, the muon spin precesses in a plane 
perpendicular to the field direction (which we define here 
as the z-direction). The time evolution of the muon spin 
polarization P{t) is determined from the /iSR "asymme- 
try" spectrum formed from the muon decay events de- 
tected in opposing positron counters 



m{t)-N2{t) 

Ni{t) + N2{t) 



aoPit) 



(1) 



where Ni(t) is the time histogram of the temporal de- 
pendence of decay positron count rate in the i*'* detector, 
and ttQ is the asymmetry maximum. In our experiments 
four positron counters were used to completely cover the 
360° solid angle in the x-y plane (see Fig.[T|). The muon 
spin polarization function for the "Left" -"Right" pair of 
detectors is defined as 



/•oo 

P^ {t) = / n{B) cos{-f^Bt 
Jo 

and for the "Up" -"Down" pair as 



(p)dB , 



(2) 



Py{t) 



n{B) cos{j^Bt 



where 7^ is the muon gyromagnetic ratio 
constant, and 

niB') = {S[B'-Bir)]), 



TT/2)dB , (3) 
is a phase 



(4) 



is the probabilty of finding a local magnetic field B in the 
z-direction at an arbitrary position r in the x-y plane. 
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FIG. 1: (Color online) Schematic of the positron counter ar- 
rangement used in the present study. The muon beam axis 
and the applied magnetic field are perpendicular to the page 
(i.e. parallel to the z-axis). The muon spin precesses in the 
x-y plane about the local z-component of the magnetic field 
before undergoing the decay /i"*" ^ e"*" -I- i/e -I- ■ In our study, 
the muon spin polarization P(t) is formed from approximately 
20 to 30 million muon decay events. The decay events de- 
tected in the "Left" and "Right" positron counters form the 
a;-component of P(t), while decay events detected in the "Up" 
and "Down" positron counters form the y-component. 



III. DATA ANALYSIS METHOD 

The TF-/iSR time spectra for each sample were fit 
assuming the following analytical solution of the GL 
equations^P for the spatial field profile of the ideal FLL 



Biv) = BoJ2- 
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where G are the reciprocal lattice vectors, Bq is the av- 
erage internal magnetic field, F{G) — uKi{u) is a cut- 
off function for the G sum, Ki{u) is a modified Bessel 
function, and u = ^/2£,abG. The cutoff function F{G) de- 
pends on the spatial profile of the superconducting order 
parameter at the center of the vortex core. Consequently, 
^(jb is a measure of the vortex core size. The FLL in all 
samples was assumed to be hexagonal. Neutron scat- 
tering experiments on fully-doped YBCO^"'^ indicate that 
the FLL below w 40 kOe is only slightly distorted 
from hexagonal symmetry due to a/b anisotropy . We 
find that accounting for this small distortion changes the 
values of S.ab and Xat by less than 5 %. Consequently, 
the FLL was assumed to be hexagonal for all samples 
studied. 

To avoid the difficulty of modelling the contribution 
of electronic magnetic moments to the ^SR line shape, 
we restricted our study to YBCO crystals free of static 
or quasistatic spins. For the applied fields considered in 
this study, this has been determined to be the case for 
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FIG. 2: (Color online) ^SR line shapes in single-crystal 
YBCO. Fourier transforms of the TF-^SR signal from y = 
6.75 and y = 6.95 samples at T = 2.5 K and = 4.92 kOe 
(green circles). The righthand peak for y — 6.95 is a back- 
ground signal coming from muons stopping outside the sample 
(Note, the background and sample peaks nearly coincide in 
the nSK line shape for y — 6.75). The red curves through the 
data are the Fourier transforms of the fits in the time domain. 
The contour plots show the corresponding spatial dependence 
of the supercurrent density j(a;, t/) = |V x B(a;, j/)|, providing a 
visual illustration of the change in core size with hole doping. 



oxygen content t/>6.50. 

To properly account for disorder of the FLL, the di- 
mensionality of the vortices must be considered. Joseph- 
son plasma resonance measurements on YBa2Cu306.5o 
ortho-II single crystals grown by Liang, Bonn and 
Hardy indicate that the vortices are 3D-like at low 
temperatures!)^ while mutual inductance measurements 
on thin films by Zuev et al show that even severely un- 
derdoped YBCO is quasi-2D only near TcJ^ Since the 
focus in the present study is on the variation of Xab and 
^(jf, in higher doped samples at low temperatures, the 



vortices are assumed to be rigid 3D lines of flux. This as- 
sumption is also consistent with the observation of highly 
asymmetric ^SR line shapes for all of our samples at low 
T (see Fig. [2]). For rigid flux lines, random displacements 
of the vortices from their positions in the ideal hexago- 
nal FLL are accounted for by convoluting the theoretical 
line shape n{B) by a Gaussian distribution of fields.^'* A 
Gaussian function also describes the local distribution of 
dipolar fields originating from static nuclear momentsi^ 
Taking into account both sources of line broadening, the 
corresponding theoretical polarization functions are 

P.W = e-'^'«*'/M n{B)cos{-ff,Bt + ^)dB, (6) 
Jo 

POO 

Py{t) = e~'"'ff*'/2 / n(S) cos(7^Si -I- - 7r/2)d£(7) 

^0 

where, 

C^cff = C^dip + <^FLL J (8) 

is an effective depolarization rate due to nuclear dipole 
moments {(Jdip) and FLL disorder (ctfll)- Values for (Jdip 
were obtained by fitting the TF-/iSR signal above Tc to 
the theoretical polarization function 

P[t) = e^'^^ip*'/^ cos(7^Si + (/)) . (9) 

To account for the background signal from muons that 
did not stop in the sample, an additional term of the form 
(1 - /)e-''B*V2 cos(7^Bi + 0b) was added to Eq. ^ and 
to Eq. ([H]), where / is the fraction of muons that stopped 
inside the sample. Values of / for the different samples 
ranged from 0.8 to 0.9. 

In the present work there are two marked improve- 
ments over the analysis done in our previous studies 
of YBCO single crystals^Si^liS^ that assumed the spa- 
tial field profile of Eq. ([S]): (i) The earlier works used 
the asymptotic limit Ki{u) — y^7T/2uexp{~u) {u ^ 1) 
for the Bessel function that appears in the cutoff F{G), 
whereas here Ki(u) was evaluated numerically, (ii) The 
second improvement is that due to increased computer 
speed, B{r) was calculated at 15,132 equally-spaced lo- 
cations in the rhombic unit cell of the hexagonal FLL, 
compared to 5,628 locations in previous works. Further 
increasing the number of real-space points sampled in the 
FLL unit cell did not result in appreciable changes in the 
fitted parameters. We note that both improvements in 
our data analysis method influence the absolute values of 
Xab and ^afc, but the temperature and magnetic field de- 
pendences of these parameters remain qualitatively sim- 
ilar to that determined in our previous studies. 

IV. RESULTS FOR Kb 
A. Temperature dependence 

Figure [3] shows the temperature dependence of Xat at 
low T determined at two different values of the applied 
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FIG. 3: (Color online) Temperature dependence of Xat at (a) 
H — 5 kOe and (b) H = 15 kOe. The solid curves are fits 
to Xab{T,H) = Xab{0,H) + aT", where a and n are field- 
dependent coefficients. 



magnetic field. The solid curves are fits to Xab{T,H) = 
Aaf,(0, H) + aT", where a and n are field-dependent co- 
efficients. The dependence of on T is shown in 
Fig. m for selected values of the applied field. An inflec- 
tion point at T w 20 K is visible in some of the lower 
field data. This feature was also apparent in our pre- 
vious measurements of YBa2Cu3 06.95-^^ Harshman et al 
have argued that the inflection point is caused by thermal 
depinning of vorticesr-^ although an invalid treatment of 
the data was used to support this assertion;^ Recently, 
Khasanov et al have ruled out depinning as the source of 
a similar inflection point in the temperature dependence 
of measured in Lai.83Sro.i7Cu04 by TF-^SRj^ In- 
stead they attribute this feature to the occurrence of 
both a large d-wave and a small s-wave superconduct- 
ing gap. As we will show later, the anomalous magnetic 
field dependence of the vortex core size in YBCO can be 
explained by an induced superconducting energy gap in 
the CuO chains that run along the b direction. Theoreti- 
cal calculations by Atkinson and Carbotte'^*' for a d-wave 
superconductor with proximity-induced superconductiv- 
ity in the CuO chains, show that 1/A^j(r) exhibits an 



inflection point caused by an upturn of 1/A^ (T) at low T 
(where Ab is the penetration depth in the b direction). 

In Fig. El it is shown that Xlb{T ^ 0) / Xl^iT) exhibits 
a near universal linear temperature dependence at low 
T. We attribute deviations from universal behavior near 
Tc to softening of the FLL, which narrows the /itSR line 
shape and enhances the fitted value of Xab- The uni- 
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FIG. 4: (Color online) Temperature dependence of l/A^j at 
H = 5 kOe. The insets show additional low temperature mea- 
surements at _ff = 2 kOe or i/ = 15 kOe. The solid curves 
through the data points are merely guides for the eye. 



versal scaling implies that A^b(T 0)TcVfZ'^^ /va is a 
constant where vf is the Fermi velocity, va is a velocity 
corresponding to the slope of the gap at the nodes, and 
Z*^ is a charge renormalization parameterizing the cou- 
pling of the quasiparticles to phase fluctuations. Using 
values of va from thermal conductivity measurements^ 
we find that Z"^ is bascially doping independent. 



B. Magnetic field dependence 

Figure [6] shows the magnetic field dependences of Xab- 
Here we stress that the observed behaviors do not im- 
ply that the magnetic penetration depth or superfluid 
density depend on field in this way. The sublinear de- 
pendence of Xab on H is primarily due to the failure of 
Eq. ^ to account for all field-dependent contributions 
to the internal magnetic field distribution. In Refsi^i^i 
the strong field dependence of Xab in YBa2Cu3 06.95 de- 
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FIG. 5: (Color online) Temperature dependence of l/\ab ^.t 
H = 5 kOe, normalized to l/A^i,(T 0) and plotted as a 
function of reduced temperature T/Tc. The solid line through 
the data is a fit given by Xlb{T 0)/\lb{T) = 1 - 0.54T/Tc. 



termined by /iSR was explained by the high anisotropy 
of the d-wave superconducting energy gap not accounted 
for in Eq. ([5]). A nonlocal supercurrent response to the 
applied field in the vicinity of the vortex cores stemming 
from the divergence of the coherence length at the gap 
nodes modifies the spatial distribution of field. With in- 
creasing H, the increased overlap of the regions around 
the vortex cores reduces the width of the /xSR line shape. 
The gap anisotropy also results in a nonlinear supercur- 
rent response to the applied field, resulting from a quasi- 
classical 'Doppler shift' of the quasiparticle energy spec- 
trum by the flow of superfluid around a vortexi^^ When 
the Doppler shift exceeds the energy gap, Cooper pairs 
are broken, and Xat increases. 

These effects are not restricted to c?-wave supercon- 
ductors. Sizeable nonlinear and/or nonlocal effects can 
also occur in s-wave superconductors with a smaller en- 
ergy gap on one of the Fermi sheets and/or a highly 
anisotropic Fermi surface. Moreover, these anisotropics 
result in a rapid delocalization of quasiparticle core states 
with increasing H that modify n(B). Indeed, strong 
field dependences of Xab fi'om Eq. ([5]) have been ob- 
served in the multi-band superconductor NbSe2r^ and 
the marginal type-II superconductor VJ^ It has been 
experimentally established for a variety of materials 
including YBCO,^ that the _ff — > extrapolated value 
of Xab agrees with the magnetic penetration depth mea- 
sured by other techniques in the Meissner phase. Conse- 
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FIG. 6: (Color online) Magnetic field dependence of the ex- 
trapolated zero-temperature value of \ab- The solid curves 
are fits to Aai,(0, iif) — Aa6(0, 0) + PVH, where the coefficient 
P decreases with increasing hole doping. 

quently, we stress that only \ab{H — > 0) can be consid- 
ered a "true" measure of the magnetic penetration depth. 



C. Hole doping dependence 

In Fig. [3 we show as a function of 1/A^j,(r -> 0) 
at two different fields. The more inclusive data set at 
iJ = 5 kOe is described by cx (1/A^f,)°-^^, which devi- 
ates substantially from the linear scaling in the Uemura 
plot. The power 0.38 is surprisingly close to 0.43 deter- 
mined by Zuev et al in a Meissner phase study of severely 
underdoped YBCO thin filmSfiS It is surprising because 
these thin films have a superfluid density that is signif- 
icantly lower than in single crystals It is known from 
microwave studies of YBCO that the doping dependences 
of Aa and A;, are not the samcy^, due to the conductivity 
of the CuO chains. '^^ Thus there is no reason to expect 
the power ~ 0.4 to be universal for the cuprates. While 
a sublinear dependence of Tc on 1/A^j, has also been in- 
ferred from more recent /xSR measurements of the muon 
depolarization rate a in other high- Tc superconductors;^ 
the contributions of magnetism and FLL disorder to the 
/iSR line shape were not factored out. 

Figure [5] shows l/A^^ as a function of hole doping, 
where the values of p are determined from the depen- 
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FIG. 7: (Color online) Dependence of Tc on l/Xlt{T 
0). The two data points on the far right are for 
the overdoped sample (Y,Ca)Ba2Cu306.98- The solid 
curves are fits to the data for YBa2Cu30j;, yielding 



Tc = (19.3 ±1.7 K-Mm^)[l/Af 



i0.38±0.02 



and Tc = (25.5 ± 



1.9 K-^im2)[l/A^i,]°-^^*° °^ a.t H = 5 kOe and H = 15 kOe, 
respectively. 
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FIG. 8: (Color online) Dependence of 1/A^6(r^0) on hole 
doping concentration p ai H = 5 kOe. The dashed curve is 
the function 1/A^i, oc (p— 0.05)^ *^, where p = 0. 5 is the critical 
hole doping concentration for the onset of superconductivity. 
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FIG. 9: (Color online) (a) Hole doping dependences of the 
Gaussian depolarization rates ctpll and crdip at T — ^ K 
and H = 5 kOe. (b) Hole doping dependence of the root- 
mean-square displacement (s^)^^^ of the vortices from their 
positions in the perfect hexagonal FLL at T — > K and H — 5 
kOe plotted as a percentage of the intervortex spacing L. 



dence of Tc on p presented in Ref.— for similar YBCO 
single crystals. The behavior is consistent with other 
studies on cuprates indicating that the maximum value 
of is not reached before p w qjq 38,39 q^j, f^^^g^ 

are in the range 0.103 < p < 0.192 and are described by 
V-^ab fx (p - 0.05)1-42, g^(,gp|- j^gg^j. p ^ Q = 1/8. 



The hole doping dependences of (Tdip and ctfll at 
H = 5 kOe are shown in Fig. IH^a). While cxdip is in- 
dependent of p, cTpLL basically tracks l/A^^. Using the 
fitted values of ctfll and Xab we have calculated the hole 
doping dependence of the root-mean-square displacement 
(5^)1/2 of the vortices from there positions in the perfect 
hexagonal FLL. As shown in Fig.jSI^b), the degree of FLL 
disorder is small and as expected highest in the Ca-doped 
sample. 
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FIG. 10: Magnetic field dependence ol the extrapolated 
zero-temperature value of £^ab- The solid curves are fits to 
(,ab{0, H) = a + h/\/H, where the coefiicients a and b depend 
on the hole doping concentration. 
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FIG. 11; Calculated field-dependent core size R relative to 
the BCS coherence length for single-layer and bilayer mod- 
els (see Appendix). Also shown is the field dependence of 
£,ab{T — »■ 0) in the y = 6.95 sample normalized to the ex- 
trapolated value ^ab{H Hc2) = 18.5 A and plotted as a 
function of H/Hc2, where Hc2 = $o/2tvS^I^{H Hc2)- 



V. RESULTS FOR 
A. Magnetic field dependence 

The field dependences of ^afc are shown in Fig.[Tni The 
increase in ^ab at low field, which corresponds to an ex- 
pansion of the vortex cores, was previously reported for 
YBasCugOg.eo (Ref.-^) and YBazCugOe.gs (Refs^^SS). 
In all samples we find that ^ah scales as 1 / ^/H, which is 
proportional to the intervortex spacing. In s-wave super- 
conductors the field dependence of ^afc has been shown to 
originate from the delocalization of bound quasiparticle 
core stateSfi^ii^ This is because a change in the spatial 
dependence of the pair potential accompanies the change 
in electronic structure of the vortex cores. In YBCO the 
low-energy quasiparticle core states should be extended 
along the nodal directions of the d-wave gap function^S 
This allows for a large transfer of low-energy quasipar- 
ticles between vortices at low field, which is further en- 
hanced by an increase in vortex density. Hence the vortex 
core size is predicted to shrink with increasing field^ii^ 
However, the field dependence of in YBCO is con- 
siderably stronger than predicted for a pure d-wave su- 
perconductor. Consequently, we consider an alternative 
explanation. 

In Fig. [n] we show that the field-dependence of £_ab, 
in particular the upturn at low field, can be explained 
by the presence of the CuO chains. The calculations of 
the core size R are based on a semiclassical Doppler-shift 
approximation (see Appendix) for either a single-layer 
model representing a superconducting Cu02 plane or a 



proximity-coupled model representing a Cu02-CuO bi- 
layer. In the bilayer model there are two distinct en- 
ergy scales for pair breaking: the energy gap associated 
with Cooper pairs in the Cu02 planes, and a smaller 
proximity-induced gap associated with the chains. It is 
the latter scale which is responsible for the expansion of 
the vortex cores at low field. 



B. Hole doping dependence 

Figure [12] shows the hole doping dependences of 
^ab{T^O) at H = 5 kOe and H=15 kOe. Qualitatively, 
the doping dependence of ^at is similar to that reported 
by Ando et al from magnetoconductance measurements 
on detwinned YBCO single crystals4^ This result is also 
shown in Fig. [121 but plotted as ^ab versus y. Note that 
our data must be plotted as ^ab versus p, because the 
hole doping concentration of our y = 6.60 single crys- 
tals (grown in a different kind of crucible than the other 
samples), is smaller than that y — 6.57 (see Table HI. 
The general trend of all data sets is an increase of ^ab 
with decreasing p. Such behavior has also been observed 
in the underdoped regime of La2-a:Sra;Cu04r^i^'^ and 
Bi2Sr2Cu06+5i'*^ With increasing magnetic field, our val- 
ues for ^ab approach those determined by Ando et al. 
Note that based on our proximity-induced model for the 
field dependence of £,ab (see Appendix) , it is the high-field 
values of ^ab that refiect the intrinsic superconductivity 
of the Cu02 planes. 

Since the doping dependences of £,ab 'At H — 5 kOe and 
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FIG. 12: Hole doping dependence ol ^abiT ^ 0) at _ff = 5 
kOe and H — 15 kOe. Also shown is data for Irom Ref.— 
plotted as a function of y, rather than p. 



H = 15 kOe in Fig. [12] are similar, we expect the hole 
doping dependence of H*2 = $o/27r^^f, to qualitatively 
resemble that of the upper critical field Hc2- Figure [T51 
shows the hole doping dependence of H*2 calculated from 
the values of ^ab at H — 5 kOe. Consistent with the data 
of Ando et al, Hc2 decreases with decreasing p in the 
underdoped regime of YBCO, and displays a dip near 
1/8 hole doping. 



VI. SUMMARY AND CONCLUSIONS 

We have simultaneously determined the hole dop- 
ing dependences of the magnetic penetration depth and 
the GL coherence length in the underdoped regime of 
YBa2Cu30y. This was achieved by fitting /uSR mea- 
surements in the vortex state to an analytical solution 
of the GL equations for the internal magnetic field dis- 
tribution. In this type of analysis the magnetic pene- 
tration depth is strictly defined as the H ^ extrapo- 
lated value of the fitted parameter Aafc. The accuracy of 
this definition was established in previous studies of con- 
ventional superconductors ] Here we have presented 
measurements showing a refinement of the Uemura plot 
for YBa2Cu30j,, where Tc is plotted as a function of the 
isolated quantity rather than the the muon de- 

polarization rate a. We find that Tc exhibits a strong 
sublinear dependence on suggesting that is not 

directly proportional to the superfluid density ps- This 
result supports the same conclusion reached in several 
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FIG. 13: Plot of H*2 = <l>o/27r^^6(r ^ 0) as a function of 
hole doping (red circles), using the values of (,ab{T — > 0) at 
H = 5 kOe. Also shown are the values of Tc (black circles). 
The black curve is the relation between T and p from Ref.—. 



recent Meissner phase studies of YBCO. 

We have reported here a reduction of 1/A^j, near 1/8 
hole doping concentration. Suppression of 1/A^^ or the 
muon depolarization rate a near 1/8 hole doping has 
only previously been observed in cuprates where p is 
controlled by cation substitution, ■^^'^^ and was believed 
to indicate a tendency toward static stripe formation. 
Static stripes over the doping range investigated here 
were recently ruled out by inelastic neutron scattering 
experiments on YBCO 4^ We ourselves find no evidence 
for static spins in zero-field /iSR or TF-/j,SR measure- 
ments on our samples. ■^^ However, the suppression of 
superconductivity near p — 1/8 could be caused by 
fluctuating stripes, recently argued to be relevant in 
YBCO and other cuprates4^ Experimental evidence for 
dynamic stripes in YBCO includes the detection of low- 
energy one-dimensional incommensurate modulations in 
YBa2Cu306.5o by inelastic neutron scatteringi^ 

Further evidence for suppression of superconductivity 
near p—1/8 is found in the hole doping dependence of 
S,ab- In our measurements ^at is a parameter that charac- 
terizes the size of the vortex cores. While it mimics the 
behavior of the GL coherence length, £^ab is large at low 
held due to the contribution of the CuO chains to the 
spatial dependence of the superconducting order param- 
eter. Enhancement of the GL coherence length or vortex 
core size near 1/8 hole doping, has also been observed 
in La2_£i;Sra;Cu044^ Calculations by Mierzejewski and 
Maska show that static or quasistatic stripes actually in- 
tensify Hc2 by reducing diamagnetic pair breakingj^ and 
hence cannot explain the growth of ^at near p=l/8. On 
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the other hand, Kadono et al have shown that an ex- 
pansion of the vortex cores with decreasing hole doping 
can result from a strengthening of antiferromagnetic cor- 
relations competing with superconductivity.^- Thus, dy- 
namic stripes are a viable explanation for the increased 
size of the vortex cores near 1/8 hole doping. 
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between sites i and j within and between layers respec- 
tively, while Ay is the superconducting order parameter 
along bonds connecting nearest neighbour sites i and j. 
From the form of Eq. (fTU]) . it is apparent that A^ only 
couples quasiparticles belonging to the plane layer. The 
single-layer Hamiltonian is obtained by setting t± — 0. 

A magnetic field H applied perpendicular to the lay- 
ers induces circulating currents in the superfluid. The 
superfluid velocity is given by Vg — • Ps where M is 
the effective mass tensor, Ps(r) — (2efi./c)A(r) -I- hW(f>{r) 
is the superfluid momentum, A(r) is the magnetic vector 
potential and 0(r) is the local phase of the order param- 
eter. In the limit that (f>{r) and A(r) are slowly varying 
functions, one can treat the superflow as uniform in the 
neighbourhood of r. Then, one can make a local gauge 
transformationSi'53 such that the phase is removed from 
the order parameter and appears instead in the hopping 
matrix elements tn^ij: 



Appendix: Semiclassical Calculation of the Vortex 
Core Size 



Calculations of the vortex core size are based on a gen- 
eralisation of the so-called "doppler-shift" approximation 
for the vortex structure^^'^^ to the case of YBCO, which 
is a multiband superconductor. In the case of YBCO, 
there is strong evidence that both the two-dimensional 
Cu02 planes and one-dimensional CuO chains supercon- 
duct. Furthermore, it is likely that the chains are intrin- 
sically normal, but are driven superconducting by the 
proximity effect. Proximity models for YBCO have been 
extensively described elsewhere^^i^. The essential idea is 
that the superconductivity originates from a pairing in- 
teraction which is confined to the two-dimensional Cu02 
planes, and that the mixing of chain and plane wavefunc- 
tions induces superconductivity in the one-dimensional 
chain layers. 

We adopt a simplified bilayer model consisting of a sin- 
gle plane and a single chain, with one Wannier orbital 
retained per unit cell for each layer. For comparison 
purposes, calculations are also performed for a single- 
layer model of an isolated superconducting plane. The 
Bogoliubov-deGennes Hamiltonian for the bilayer is 



f ■ — f . .p-*Ps(r)-(ri-rj)/2?i 



H 



tl,ij Aij 


















t^5^,j 









-t±^i,j 





~t2,ij* . 





where = [■«/'LT'^iii^2iT^2ii] and ^p^^^ (V'mi) are cre- 
ation operators for quasiparticles (quasiholes) at lattice 
site i in layer n. Here, we take n = 1 for the plane 
layer and n = 2 for the chain layer. The parameters tnjj 
and tj_5i_j are the single-electron hopping matrix elments 



tn,i] + -z^ij ■ Ps(r) 



(11) 

(12) 



where t^^ij are the hopping matrix elements in zero-field 
and Vij = —itn,ij{Yi — rj)/h are the matrix elements of 
the zero-field quasiparticle velocity. Equation fol- 
lows from Eq. I|lip in the limit that Ps is small. Then, 
the order parameter takes on the simple d-wave form 
Ay = iA(— 1)^'~^J which, in reciprocal space, corre- 
sponds to Ak ~ A[cos(/ca;a) — cos(fcya)], where a is the 
lattice constant. The local gauge transformation leads 
to a doppler-shifted spectrum and is exact in the limit 
of slowly varying superfluid velocity. This procedure has 
been shown, in many circumstances, to provide a reason- 
able description of the vortex lattice^iS. 

We take band structures which are appropriate for 
YBCO and adopt 



tlSj — 



to i^j 

ti i,j are nearest neighbours 

^2 hj are next-nearest neighbours 



(13) 



and 



t2.i] — 



'3 . . . (14) 

^4 i, j are nearest neighbours along y 



For this work, we measure energies in units of \ti\ and 
take {tQ,...,t4} = {1,-1,0.45,2,-4}. In reciprocal 
space, the dispersions of the isolated plane and chain 
layers are then eik = to + 2ti[cos{kxa) + cos(fcj,a)] -I- 
4^2 cos(fc^a) cos(fcj^a) and e2k = ^3 + 2^4 cos(fcj,a), respec- 
tively. The chain-plane hopping matrix element t± is not 
well known in YBCO and is taken to be t± = 0.75. 
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For a slowly varying Ps (r) we can locally Fourier transform the Hamiltonian in the neighbourhood of r to give 



k 



ivife-ps(r) Afc(r) 
Afe(r) -eifc + ^vife • ps(r) 



-tj_ 







(2k 



5V2fc • Ps(r) 












-f2k 



\^2k ■ Ps(r) 



(15) 



where v„k — h^^denk/dk and ~ N^^/^ J2i where 
N is the number of fc-points in the sum in Eq. (jisp . 

We need to make an ansatz for Ps(r)- For a single vor- 
tex in an isotropic medium, one has Ps(r) — {2nh/r)9, 
where 9 is the azimuthal unit vector and the radius r 
is measured relative to the centre of the vortes*^. For 
the bilayer model, however, Ps(r) is not isotropic: the 
chains provide a conduction channel along the y direc- 
tion which is in parallel with the isotropic plane con- 
duction channel. We mimic this anisotropy by assum- 
ing that the superfluid momentum will be similar to 
that of a single-layer superconductor with an anisotropic 
(diagonal) effective mass tensor M with Alyy < M^x- 
(For the single-layer model, we take Mj.x = ^yy ) We 
then have two requirements which must be satisfied: 
Vxps = 2i:hY,Ti^'^{-r-^) and V-v^ ee V-M^i-ps = 0. 
The first requirement introduces vortex cores at the vor- 
tex lattice sites R, while the latter incompressibility re- 
quirement is strictly true in regions where A(r) is uni- 
form. This pair of equations is solved by 



im-{Gx z) 



G 



(16) 



where indicates that G = is excluded from the 
sum, G are reciprocal lattice vectors of the magnetic unit 
cell (we assume a square lattice here) with area and 
magnetic length L. The results do not depend strongly 
on the ratio Myy/Mxx, which we take to be 0.6 for the 
parameters chosen above. This choice minimizes V • j (r) , 
where j(r) is the total (plane and chain) current in the 
bilayer. 



j(r) 



1 ^ 



Vlk 



V2k)r 



(17) 



and (. . .)r indicates the expectation value with respect 
to H{r), Eq. p^ . In principle, one could improve on 
the approximation of Eq. (fT6|) by determining Ps(r) self- 
consistently from j (r) ; however this will not change the 
qualitative physics of the vortex core expansion. 

We then solve self-consistently for the order parameter 



V 

A(r) = 2_^[cos{kxa)-cos{kya)]{^/Jl^kl'^Plk^)r, (18) 



with V — 1.7. Self-consistent solutions find that A(r) 
vanishes near the vortex core centre and obtains an 
asymptotic value Amax = 0.35 far from the vortex core. 
In order to measure the vortex core size, we define a 
quantity (5A(r) = Amax — A(r). The vortex core size is 
then defined by the first moment of the radial coordinate 
r with respect to (5A(r): 



Er'^'^A(r) 
Er'5A(r) ■ 



(19) 



where r = corresponds to the vortex core centre. For 
presentation purposes, R is shown relative to the BCS co- 
herence length ^0 = fivp /T^^max, where vp the average 
of the Fermi velocity on the Fermi surface. The magnetic 
field is related to the magnetic length by H = ^q/L^ 
where $o is the superconducting flux quantum. For pre- 
sentation purposes, H is scaled by the upper critical field, 
H,2 = ^o/2t:^1 so that H/H,2 = ^ttQ/L^ 
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